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Abstract 



o 

» ! ' Vortons can be viewed as (flat space-) field theory analogs of black rings in general relativity. They 

are made from loops of vortices, being sustained against collapse by the centrifugal force. In this work 
we discuss such configurations in the global version of Witten's (7(1) x (7(1) theory. We first consider 
solutions in a flat spacetime background and show their nonuniqueness. The inclusion of gravity leads 
to new features. In particular, an ergoregion can occur. Also, similar to boson stars, we show that 
the vortons exist only in a limited frequency range. The coupling to gravity gives rise to a spiral-like 
frequency dependence of the mass and charge. New solutions of the model describing 'semitopological 
, vortons' and 'di-vortons' are also discussed. 

cr 



^jrj! 1 Introduction 



Vortons are certainly among the most interesting particle-like solutions in field theory. Within a rather 
informal definition, they are ring-like solitons in some nonlinear field theory which is supposed to possess 
(Zf^ , also string-like vortex solutions. Then the heuristic construction of a vorton involves taking a piece of finite 

length of such a vortex, bending it and gluing its ends to form a loop. The resulting configuration -the 
vorton, necessarily caries a nonzero angular momentum, such that the tension of the loop is compensated 
by the centrifugal force. This stationary ring does not radiate classically, and at large distances looks like a 
CO '. point particle with quantized charge and angular momentum. 

Vortons were initially proposed 25 years ago [1, 2, 3] in the context of the local (7(1) x (7(1) theory of 
superconducting cosmic strings of Witten [4] . This model possesses vortex solutions with a non-dissipativc 
current flowing along the string due to a complex scalar field which vanishes asymptotically. Vortons are 
closed loops of such superconducting strings, carrying both current and charge. However, vortons may exist 
in a variety of other theories as well, ranging from cosmology [5] to QCD [6] and condensed matter systems 
[7]. One can speculate that vortons may perhaps even exist in the Standard Model of Particle Physics. 

The study of vortons has been a fruitful avenue of research since its inception. Some of the ideas which 
appeared in this context have found interesting applications in other domains of physics. For example, 
vortons share the same heuristic construction with black rings in d > 5 general relativity [8, 9], with vortices 
playing the role of black strings there. Moreover, a general formalism proposed by Carter in [10] to describe 
the properties of large radius vortons has found recently beautiful applications and generalizations in the 
physics of blackfolds in higher dimensional Einstein gravity [11]. 

At the same time, different from black rings in general relativity 1 , there are no closed form vorton 
solutions. Moreover, even the numerical study of such configurations is difficult, due to the existence of 
multiple scales in the problem and the complicated behaviour of the matter functions. 

In fact, to our knowledge there is no explicit construction in the literature of vorton solutions in a model 
containing gauge fields. However, by now, there are a number of independent studies of vortons in the global 
version of Witten's (7(1) x (7(1) model, which contains only two complex scalars. The numerical study of 



1 Note, however, that the black rings are known in closed form for d = 5 spacetime dimensions only [8]. The d > 5 black 
rings were studied based on analytical solutions valid in the large radius limit [9] or on numerical nonperturbative solutions 
[12]. 
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such solutions has started with the work of Lemperier and Shellard in Ref. [13]. There, the full hyperbolic 
evolution problem for the fields, with the initial data corresponding to a vortex loop has been considered. 
The results in [13] suggested the existence of vortons which are stable against perturbations. However, not so 
many details on the properties of such configurations were given there. In addition, Lemperier and Shellard 
did not consider precisely the rigid version of Witten's model, but, in order to improve the numerics, added 
an extra Q-ball type interaction term. 

The existence of stationary, axially symmetric vortons in Witten's model has been reported in Ref. [14]. 
There, the solutions have been constructed by solving a set of elliptic partial differential equations with 
suitable boundary conditions 2 . The configurations in [14] possess a nonzero value of the angular momentum 
and have ring profiles, describing thick vortons. Also, the results there point towards the existence of 
some analogies between vortons and spinning Q-balls in a simpler model with a single scalar field with a 
non-renormalizable potential [16, 17]. 

Large radius vortons have been reported in [18] for a different choice of the parameters of Witten's 
model and a different numerical approach than the one employed in [14]. Ref. [18] has investigated also the 
numerical stability of such configurations with respect to non-axially symmetric perturbations. The results 
there show that (at least) thin vortons possess a pinching instability. 

The main purpose of this work is to enlarge our understanding of vortons by including the effects of 
gravity, at a nonperturbative level. This is a legitimate task since the geometry in the vorton's core may 
exhibit large deviations from the flat spacetime limit. Also, as shown by a number of other models [19, 20, 21] , 
the inclusion of gravity usually leads to a variety of new effects. Indeed, as we shall see, this is also the case 
for vortons. For example, they may possess an ergoregion, with the associated instability. Also, gravitating 
vortons exhibit a rather different type of frequency dependence. 

However, as a first step towards this task, we shall start with a study of vortons in a flat spacetime 
background. The solutions are easier to study in this limit, while some of their basic properties remain 
unaffected by the inclusion of gravity. Based on the heuristic construction of vortons as pieces of vortices, we 
propose a classification of the resulting objects in terms of an integer k related to the winding along a torus 
generator. Then it follows that vortons can be viewed as the second member (k = 1) of a general family 
of solutions. An unexpected result reported here is the nonuniqueness of vortons, that is the existence of 
different solutions for the same input parameters. Apart from solutions already known in the literature, new 
types of configurations, in particular those describing di- vortons, are also shown here. 

This paper is organized as follows. In Section 2 we recall Witten's gravitating global model, and give a 
precise formulation of the problem we propose to solve. We continue in Section 3 with a discussion of the 
properties of vortons in a flat spacetime background. The effects of gravity on the properties of vortons are 
discussed in Section 4. Two kinds of solutions are described there, corresponding to 'semitopological' vortons 
and 'true' vortons. Section 5 gives our conclusions together with a discussion of possible generalizations and 
open problems. In the Appendices we present the system of differential equations solved in this work (A), 
the components of the energy-momentum tensor (B), and the spherically symmetric limit of our solutions 
(C). 

2 The model 

2.1 The action and field equations 

Witten's ungauged model of superconducting cosmic strings [4] contains two complex scalar fields <j> and a, 
with the following Lagrangian density 3 

L s = 0^*0^ + d^d^a + 1/(101, |a|), (2.1) 

2 An independent construction of these solutions has been reported in [15], based on a very different treatment of the 
numerical problem. 

3 Note that, as opposed to previous studies [13, 14, 18], we are using in this work a mostly plus signature for the spacetime 
metric. 
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Figure 1. The potential t/(|0|, |er|) of the model is shown for a set of parameters with = rj a = 1, A^, = 41, 
A CT =40, 7 = 22. One can notice the existence of a single global minimum at \<f>\ = rjj,, \a\ = 0. 



where J7(|0|, |er|) is the scalar fields potential 



um, \a\) = \\m 2 - vD 2 + \xm 2 (w\ 2 - H) + 7i0i 2 ki 2 . 



(2.2) 



In the above potential, A^ , A CT , 77^, r] a and 7 are positive constants which are not specified apriori. A 
minimal value of the potential is achieved for \tp\ — rj^, \a\ — 0, in which case i/(|<£|, |er|) = (see Figure 1). 
This minimum is globally unique and stable if 



4 7 2 > A^A^ and 777^ > i \ a rf a . 



(2.3) 



Then the perturbative spectrum of the field excitations around the vacuum state \<f>\ = 77^, \a\ = consists 
of two massless Goldstone particles, corresponding to excitations of the phases of the fields, and of two Higgs 
bosons with the masses 



Ma, 
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When coupling to gravity, the dynamics of the solutions is governed by the following action 



d 4 x y 



( R 



L, 



(2.4) 



(2.5) 



where g is determinant of the spacetime metric <? M „ , G is Newton's constant and R is the curvature scalar. As 
usual in general relativity, variation of the action with respect to the metric leads to the Einstein equations 



E^v — R^ v — -g^R — 8itG T^ lv — , 



(2.6) 



with the stress-energy tensor 

dL s 



g^vL s 



= (4,*^ „ + 4>* v 4 iU , + a^a tV + a* u a „) - 9fil/ [d a (t>*d a 4> + d a a*d a a + U{\4>\, \a\j\ . (2 
Variation of (2.5) with respect to the scalar fields leads to the matter equations, 



(2.7) 



□ 



dU 







□ ■ 



dU 



(7 = 



(2.9) 
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where □ represents the covariant d'Alembert operator. 

We close this part by noting the existence of several interesting limits of the considered model. First, for 
a real field <f> and a choice X a = in the potential (2.2), one recovers the model proposed by Friedberg, Lee 
and Sirlin (FLS hereafter) in [22], with 



+ d^d^a + -A (H 2 - 4? + 7 |0| 2 M 2 . (2.10) 



As shown in [22], this model possesses spherically symmetric, nontopological soliton solutions. In Section 
3 we shall give numerical evidence for the existence of spinning, axially symmetric generalizations of those 
configurations. 

Another important limit is found in the special case r\ a — r/^ by redefining the parameters in the potential 
as X,p = \ a = (3, and f] — \fi ^ 7o- Then, after taking the limit /? — > oo, the theory (2.1) becomes a sigma 
model for two complex fields <f), a with 



L s =d„<j>*d»c]) + d„a*d»o + 7oM>n (2-11) 



subject to the constraint 



H 2 + M 2 ^. 

Note that this model has less degrees of freedom and no free parameters, since the value of 7 in (2.11) can 
be changed by rescaling the spacetime coordinates. Vorton solutions of this model (although for a different 
formulation) have been discussed in [7]. 

Finally, one should remark that the system described by (2.5) possesses an interesting limit with cf) = 77^, 
i.e. one scalar field only: 

L s = 8^*8^0 + Ml\a\ 2 + ^X a \a\ 4 . (2.12) 

Different from the above cases, no flat spacetime solitons exist in this case 4 . However, physically interesting 
solutions are found when gravity is included. These are the gravitating boson stars with a quartic self- 
interaction potential originally discussed in [23]. The field a then necessarily possesses a harmonic time 
dependence (note that the limit \ a = is allowed for such solutions). Their existence can be attributed 
to the balance between the dispersive effect due to the wave character of the complex scalar field and the 
self- gravitating effect. 

2.2 Global charges 

Our solutions are globally regular (i.e. without an event horizon or conical singularities), and stationary. 
Also, they approach asymptotically the Minkowski spacetime background. Then their mass M and angular 
momentum J can be obtained from the respective Komar expressions [24], 



M 



= 4^ J R^n^W, ^ = -g^/ vY<iV. (2-13) 

Here E denotes an asymptotically flat spacelike hypersurface, is normal to E with n^n^ = —1, dV 
is the natural volume element on E, £ denotes an asymptotically timelike Killing vector field and r\ an 
asymptotically spacelike Killing vector field (corresponding to azimuthal symmetry) [24] . After replacing in 
(2.13) the Ricci tensor by the stress-energy tensor (via the Einstein equations (2.6)), one finds 

M = 2 £ (t^ - l - n^"dV , (2.14) 

J = - J (t„ v - l - 9tlv Tj) n^rfdV . (2.15) 



4 Here we consider strictly positive values for the parameter A CT . However, flat spacetime solitons are found for a potential 
unbounded from below, \ a < 0, in which case they are likely to be unstable. 
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The Lagrangian (2.1) of the scalar fields has a global U(l) x U(l) symmetry, — > e tai (j), a — > e m2 cr, 
which leads to two conserved conserved currents 

j^=-i(^4>- <j>d»<j>*) , jfa=-i{v*d»v -*&>**) , with ^=0,^=0. (2.16) 

This implies two conserved Noether charges, 

Qw=fj^dV, Q {a) = JjfanpdV. (2.17) 
2.3 The axially symmetric ansatz 

The line-element for the solutions in this work possesses two commuting Killing vector fields, £ and r], with 

£ = dt , v = d v , (2.18) 

in a system of adapted coordinates. General relativity solutions with these symmetries are usually studied 
within a metric ansatz 

ds 2^ e 2U(p,z) (^k{p,z)^ df) 2 +dz 2^ + p2^^ d ^ _ e - 2 ^.*)((ft + f2(p,«)d ¥ j) 2 , (2.19) 

where (p, z) correspond asymptotically to the usual cylindrical coordinates. In the (electro-)vacuum case, it 
is always possible to set P = p, such that only three independent metric functions appear in the equations, 
and (p, z) become the canonical Weyl coordinates [25]. For more general matter sources, however, the generic 
metric ansatz (2.19) with four independent functions is needed. A variety of interesting general relativity 
solutions have been constructed within this framework [25], some of them describing toroidal configurations 5 
(see e.g. [27, 28]). 

One might wonder whether in the numerical study of gravitating ring-shaped solitons other coordinate 
systems would be better suited. For example, the use of ring coordinates seems to be an obvious choice, 
since the required symmetry of the configurations is intrinsic in this case. However, even in the absence of 
back reaction, we have met severe numerical difficulties in the implementation of a numerical scheme for 
these coordinates 6 . 

The solutions in this work are constructed by using 'quasi-isotropic' spherical coordinates (r,0), defined 
by the coordinate transformation in (2.19) 

p = rsin8, z~rcos8 7 (2.20) 

and possessing the usual range < r < oo, < 8 < it. The corresponding line-element 7 has an ADM-like 
form and reads 

ds 2 = ^(dr 2 + r 2 d9 2 ) + -v 2 sin 2 9(dip - —dt) 2 - fdt 2 . (2.21) 
/ / r 

The four metric functions f, I, m and W are functions of the variables r and 8 only, chosen such that 
the trivial angular and radial dependence of the line element is already factorized. The relation between 
the metric functions in the above line-element and those which enter (2.19) is / = e 2U P 2 / (e 4U P 2 — fi 2 ), 
Ism 2 8 = P 2 , m = P 2 /(e iu P 2 - fl 2 ) and W/r = n/(e iU P 2 - fl 2 ). The symmetry axis of the spacetime 

J Note, however, that the known toroidal solutions in four dimensional general relativity do not include black holes. Indeed, 
for d = 4 there are a number of no go results forbidding, under very general assumptions, a nonspherical topology of the horizon 
for an asymptotically flat black hole [26]. 

6 For example, the whole of spatial infinity maps in this case to a single point, which leads to bad numerical convergence. 
Also, the ring radius is an input parameter which should be specified apriori; however, this did not appear to be possible within 
our treatment of the problem. 

7 A similar line element (in particular the same metric gauge choice for the [r, 0)-part of the metric) has also been used in 
previous numerical work on spinning solutions in general relativity, starting with the early work [29, 30]. 
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is given by rj = 0. It corresponds to the z-axis with 9 = 0, w. The Minkowski spacetime background is 
approached for r — > oo, with f = I = m = 1, W = 0. 

The ansatz for the scalar fields which we shall use is similar to that proposed in the previous work [14], 
with 8 

4>= X(r,0) + iY(r,0), a = Z{r,9)e^ nv+Wt ^ (2.22) 

where X, Y, Z are real functions, and w and n are real constants. Single- valucdness of the scalar field a 
requires 

a(tp) = er(27r + tp) , (2.23) 

thus the constant n must be an integer, i.e., n = 0, ±1, ±2, ... . As usual in the literature, we shall refer to 
n as rotational quantum number. Without any loss of generality, we shall take in this work w > 0, n > 0. 

We close this part by noting that £ = d t and r\ = d v are not symmetries of the solutions in this work, 
the system being only effectively stationary 9 . 

2.4 The equations and the issue of the constraints 

After inserting the ansatz (2.21), (2.22) into the field equations (2.6), (2.9) we find a system of seven coupled 
partial differential equations that needs to be solved. These are four equations for the metric functions 
/, I, m, W (which are found by taking a suitable combination of the Einstein equations E^ + E® = 0, E£ = 0, 
El = and E l = 0, see Appendix A) and three equations for the functions X, Y and Z describing the scalar 
fields. The equations are of the form 

V 2 Ti = Ji (2.24) 

where T = (f,l,m,W;X,Y,Z), V 2 = -p-d r (r 2 d r ) + 2 . e de{sm.6de), and Ji are 'source' terms depending 
on the functions T and their first derivatives (the explicit form of the equations is presented in Appendix 
A). 

Apart from these, there are two more Einstein equations Eg = 0, E^ — Eg = 0, which are not solved 
in practice. Following an argument originally proposed in [35], one can, however, show that the identities 
\7 fl E fJ-r = and V^E^ 6 = 0, imply the Cauchy-Riemann relations 

d f V 2 + dgV! = 0, dfVi - d e V 2 - 0, (2.25) 

with Vi = V^gE^, V 2 = ^f=gr{E r r -El)l2 and dr = ^. Therefore the weighted constraints E r g and E r r -E e e 
still satisfy Laplace equations in (f , 0) variables. Then they are fulfilled, when one of them is satisfied on the 
boundary and the other at a single point [35]. From the boundary conditions (2.34)-(2.38) we are imposing, 
it turns out that this is the case for our solutions, i.e. the numerical scheme is consistent. 

2.5 Mass, angular momentum and charge 

The mass M and the angular momentum J can be evaluated based on the general Komar expressions 
(2.13). There we insert the metric ansatz (2.21), with unit vector n = —JJdt, and volume element dV = 

8 Note that (2.22) is not the general ansatz compatible with the symmetries of the line element (2.21). In principle, one can 
add an extra (r, #)-dcpcndent phase to the scalar field a by taking 

a= |<T(r,6»)|e i (*( r ' <, )+ rl ^+ 1 "') = {Z 1 (r, 9) + iZ 2 (r, 0)) e i(»¥>+™*) 

with Z\, Z2 real functions. (Solutions with such a complex scalar field were constructed in [14, 34] for a Q-ball model in a fixed 
Minkowski spacetime background.) Also, one can easily verify that Z2 = is a consistent truncation of the full model. 

9 Thus, one can say that the configurations in this work do not possess Killing vectors. The scalar field a depends in a 
nontrivial way on the Killing coordinates ip and t, though the corresponding energy-momentum tensor is still compatible with 
the spacetime symmetries. This is, however, a generic feature of d = 4 gravitating spinning configurations possessing Q-ball 
features (see also Refs. [31, 32] for examples of d = 5 boson stars and black holes with only one Killing vector). 
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V V7 V~9 dO dip. Also, we make use of the fact that both R\ and R can be expressed as total derivatives: 
V~aR\ = i- (r 2 sinfl^ (sin 2 0lW(d r W -—)- fd r f) ) + ^ ( sin^ (sin 2 9lWdgW - fdgf) J , 



dr \ 2/ 2 V r '■''■') 86 \ 2/ 2 

_. d /r 2 ; 3 / 2 sin 3 tf Tir \ 9 /rZ 3 / 2 sin 3 fl T \ 

5< = «Z ^ (r^-T^) +- oJ , 2 cW , (2-26) 



<9r V 2/ 2 v 'J 89\ 2/ 

(where -y/^g = r 2 sintf-^y 2 ). Then it follows that M and J can be read off the asymptotic expansion of the 
metric functions / and W, respectively 

/=l-^ + (l), W-™ + 0£). (2.7) 
i.e. after putting (2.26) into (2.13) one finds 

M = -5- lim r 2 d r f , J = ^- lim r 2 W^ . (2.28) 

2G r-s-oo 2G i — >o° 

Alternatively, the mass M and the angular momentum J can be obtained by direct integration of the 
expressions (2.14) and (2.15), 

M = J (2T; - S; 17) n„ - I (2 T* - T£) ^=5 dr d# , (2.29) 

corresponding to the Tolman mass, and 

J=- j T$,y/=gdrdBdtp. (2.30) 

The explicit form of the nonvanishing components of T£ for the ansatz (2.21), (2.22) is given in Appendix 
B. 

As discussed above, the global U(l) x U(l) symmetry of the scalar fields may lead to further nonvanishing 
global charges. For the ansatz of the scalar fields (2.22) with the scalar field <j) possessing no time dependence, 
there is only one nonvanishing conserved scalar charge 

Q(a) = - J j\ a )V-9drd6dip 

= 4n J™ dr £ d8 r 2 sin8^(w + n^j Z 2 = Q . (2.31) 

However, one can easily verify that the angular momentum density T* and the temporal component of the 
current jf> are proportional, Tj£ = njty Then, from Eq. (2.30) for the angular momentum J and Eq. (2.31) 
for the scalar charge Q, one obtains the important quantization relation for the angular momentum, 

J = nQ. (2.32) 

Thus Q is not an independent global charge, and the vortons' angular momentum is carried by the scalar 
field a only. Also, it follows that all solutions with J 7^ 0, thus in particular the vorton solutions, cannot be 
viewed as perturbations around some static configurations 10 . (Note, that (2.32) holds for solutions in a flat 
spacetime background as well.) 

We close this part by noting that gravitating vortons have no horizon. Therefore they are zero entropy 
objects, without an intrinsic temperature. The first law of thermodynamics [33] reads in this case 

w 

dM = wdQ = -dJ . (2.33) 
n 



10 A similar result has been found for Q-balls and boson star solutions [29, 17, 30]. 
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2.6 Boundary conditions 



In our approach, the solutions are constructed by solving the set (A.1)-(A.6) of seven elliptic partial dif- 
ferential equations supplemented with suitable boundary conditions on the boundaries of the domain of 
integration (see the discussion in Section 2.8). The choice of appropriate boundary conditions must guar- 
antee that the solutions are globally regular and asymptotically flat, and that they possesses finite global 
charges with finite densities 11 . 

Thus, at the origin r = we require 

<9r/|r=o = d r l\ r=0 = d r m\ r=Q = W\ r = = , d r X\ r=0 = Y\ r=0 = Z\ r=Q = 0, (2.34) 

The boundary conditions at infinity are such that the ground state of the model is approached, 

f\r— »oc = l\r^>oo = ?Tl|r— >oo = 1? I^lr^oo = ; -^|r— »oc = V(f>i ^Ir^oo = ^|r— >oo = 0- (2.35) 

On the symmetry axis (9 — and 9 = ir), we require the boundary conditions 

def\e=o,TT = d g l\0= o ,TT = d g m\ g=0 ^ = d W\e= o ^ = 0, d g X\ g=0 ^ = d g Y\ g=0 .^ = Z\ 0=Oi7T = 0. (2.36) 

The elementary flatness condition [25] (i.e. the absence of conical singularities) imposes an extra condition 
on the symmetry axis 

l\d=o,n = m\e=o,T!- (2.37) 

The solutions in this work possess an extra symmetry w.r.t. reflection in the equatorial plane z = 
(9 = tt/2). Namely, all functions except Y are invariant under the transformation 9 — > n — 9, while Y 
changes the sign. We use this symmetry to reduce the domain of integration to [0, oo) x [0,7r/2]. Then the 
boundary conditions which are imposed in the equatorial plane are 

def\e=iv/2 = d g l\ g=n/2 = d g m\ g=7r/2 = d g W\ g=7r/2 = 0, d g X\ g=7r/2 = Y\ g=7r/2 = d g Z\ g=v/2 = 0. (2.38) 



2.7 Scaling properties and input parameters 

A simple inspection of the field equations reveals that the model possesses seven input parameters: 
{A^, X a , 7^0, ric, 7} and {w,n}. However, as usual, some of them can be absorbed into a redefinition of 
variables together with suitable rescalings. 

Let us start with the simplest case of nongravitating solutions. Then the scaling of the fields 

<f> -> T700 , a -> 770(7 (2.39) 

together with the redefinitions 

A*->A^/»#, X a ^\ah% 7^7/^0, (2-40) 

show that one can set 770 = 1 in all field equations without loss of generality. Then 77^ appears as an overall 
factor in front of the (rescaled-) Lagrangian density. 

In the next step, we note that a rescaling of the spacetime coordinates 

jc" -> x»/L (2.41) 

(with L > an arbitrary constant) together with a redefinition of several parameters in the potential 

A0^A L 2 , X a ^X a L 2 , j^-fL 2 , (2.42) 

1 1 These boundary conditions are also implied by a power series expansion of the solutions on the boundaries of the domain 
of integration. 
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leaves the field equations invariant. This symmetry is used in this work to set 12 the length scale of the 
vortons as given by the mass of the scalar field <fi, 



L = 



Ms 



(2.43) 



This directly reveals the existence of three essential, dimensionless parameters of the nongravitating vorton 
model 



o _ R _ X ° a _ 7 

Pi - , P2 — T— , P3 - T— 

V<p A A 

From (2.3), the parameters ft should satisfy the conditions 



(2.44) 



4/3 2 >ft, /3 3 > -ftft 2 . 



(2.45) 



Note that the scalings above are independent of any ansatz. However, when considering a scalar a of the 
form (2.22), one has to take also 



w 



(2.46) 



such that the phase of a remains invariant under the scaling (2.41), with L given by (2.43). 

The scalings discussed above hold also when gravity is included. However, in that case, the presence of 
Newton's constant G leads to the occurrence of one more dimensionless coupling constant, which, following 
previous work on solitons coupled to gravity [21], we shall call a, with 



a 2 = WttGt] 2 . 



(2.47) 



Thus the constant a corresponds to the ratio between the symmetry breaking scale and the Planck mass 
Mpi = 1/VG. On general grounds, we expect this ratio to be very small; however, to understand the general 
pattern of the solutions, we shall attempt to investigate the full range of a. 

Then, for the specific ansatz (2.21), (2.22) with the above scalings, the equations for gravitating vortons 
(which are effectively solved) can be derived by extremizing the following reduced action 



Sred = J drd9 (Cg - a 2 C s ) , 



(2.48) 



with (here we define (VA) • (VB) = A. r B. r + \A fi B fi ): 



r 2 sin 



^-(V0-(V m)+ , 2 sin^ 
1 „ cot 6», , 1 



-j(l,r + 

rl r 



■l,e) H ("V 



(vwo : 

cot6» 



^(v.f) 2 



rm 



m fi) 



(2.49) 



r sm( 



f 



■^-((VI) 2 + (VY) 2 + {VZ) 2 ) + \{X 2 + Y 2 
m 4 



1)2 + ?1Z 2 (Z 2 



2p\ 2 ) 



+ /3 3 Z 2 (X 2 + Y 2 ) -hw+ ^) 2 Z 2 + 
J r 



r 2 sin 2 9 I 



Lz 2 



In the above relations we use a dimensionless radial variable r and a dimensionless frequency w given in 
units set by M^. These conventions we shall use for all numerical data reported in this work. Also, in all 
plots, the total mass M is given in units of A-kM^/X^, while the total angular momentum J and the charge 
Q are given in units of An/X^. The components T\ and of the energy momentum tensor are given in 
units of Mi/Xff, and Mj/A^, respectively. 



12 Note that other choices are possible. For example, one can also set L 2 = l/A CT r;^ or L 2 
(2.43) leads to a length scale with a transparent physical meaning. 



I/777?. However, the choice 
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2.8 The numerical method 

The set of seven coupled non-linear elliptic partial differential equations for the functions T = 
(X,Y, Z; f,l,m,W) has been solved numerically subject to the boundary conditions (2.34)-(2.38). In a 
first stage, a new compactificd radial variable x is introduced instead of r, such that the semi-infinite region 
[0,oo) is mapped to the finite interval [0,1]. (This avoids the use of a cutoff radius.) There are various 
possibilities for such a mapping, a flexible enough choice being x = r/(l + r). This involves the following 
substitutions in the differential equations 

rT, r — > (1 — x)T tX and r 2 J\ rr — > (1 — x) 2 T iXX — 2(1 — x)T tX . (2.50) 

Next, the equations for these functions are discretized on a grid in x and 0. Here we have considered various 
grid choices, the number of grid points ranging between 280 x 30 and 150 x 70. The grid covers the integration 
region 0<x<land0<6>< tt/2. 

All numerical calculations have been performed by using the professional package FIDISOL/CADSOL 
[36]. This code requests the system of nonlinear partial differential equations in the form 
P(x,6\T\T x ,Te\T x e,T xx ,TQe) — 0, subject to a set of boundary conditions on a rectangular domain. 
Besides that, FIDISOL/CADSOL requires the Jacobian matrices for the equations w.r.t. the functions T 
and their first and second derivatives, the boundary conditions, as well as some initial guess for the functions 
J ' . The solver uses a Newton-Raphson method, which requires to have a good first guess in order to start a 
successful iteration procedure. This software package provides also error estimates for each function, which 
allows to judge the quality of the computed solution. The typical numerical error 13 for the solutions reported 
in this work is estimated to be of the order of 10 -3 . A detailed description of the numerical method and 
explicit examples arc provided in [36] (see also e.g. the Appendix in [37] for a discussion within a physical 
problem). As a further check of numerics, we have verified that the families of solutions with a varying 
frequency satisfy with a very good accuracy the first law of thermodynamics (2.33). 

In this scheme, there are six input parameters: the parameters /3j in the potential of the scalar fields, the 
(scaled) frequency w and the winding number n in the ansatz for the scalar field a, and also the constant 
a as given by (2.47). The number of nodes k of \4>\ in the equatorial plane (which provides a classification 
of the solutions, see the discussion below), as well as the global charges (mass and angular momentum) are 
computed from the numerical solution. 

In constructing gravitating solutions, we start with flat spacetime configurations as initial guess (i.e. 
taking f = l = m=l, W = 0). Then we slowly increase the value of a. When the iterations converge 
we obtain, by repeating this procedure, solutions for higher values of a. A similar procedure is used to 
investigate the dependence on other parameters of the problem, starting instead with solutions with a given 
a. Also, in some of the calculations, we interpolate the resulting configurations on points between the chosen 
grid points, and then use these for a new guess on a finer grid. 

As a general remark, we have found the study of vortons more difficult than a number of other problems 
which we studied with a similar approach. First, the present system possesses a large parameter space, and 
the solutions exist only in some regions there (and thus not for arbitrary values of {(3i, w, n, a}). Moreover, 
rather good initial profiles are usually required for the convergence of the numerical iteration. Next, the 
profiles of the scalar functions exhibit large gradients in a small region around the equatorial plane, which 
requires a careful choice of the grid. Finally, the problem is further complicated by the possible existence of 
several different solutions for the same input parameters (see the discussion below). 

3 Solutions in flat space: nongravitating vortons 

3.1 The heuristic construction of a vorton and general features 

As mentioned in the introduction, vortons are essentially loops made of pieces of vortices. Remarkably, this 
simple heuristic construction leads to a number of predictions about their qualitative properties. This is why 

13 However, the errors increase dramatically when studying solutions close to the limits of the domain of existence in the 
parameter space. 
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Figure 2. The heuristic construction of a vorton starting with a piece of a vortex. 

we shall start with a discussion of this point (more details on that, including relevant plots, are given in 
Section 6 of Ref. [14]). 

As discussed by several authors (see e.g. [2, 38, 39]), Witten's model of superconducting strings admits 
cylindrically symmetric, vortex-type solutions. For a cylindrically symmetric line element with ds 2 = dp 2 + 
p 2 di/j 2 + dz 2 — dt 2 , the appropriate ansatz for the scalar fields reads 

0=|0(p)|e^, a=\a{p)\e i ^ z - wt \ (3.51) 

As usual, the single- valuedness of the scalar fields imposes k to be an integer, k = 0, ±1, . . . , while there are 
no apriori restrictions for k z ,w in (3.51). 

A study of the asymptotic form of the solutions reveals that the scalar field 4> approaches the vacuum 
expectation value (v.e.v.) for p — > oo; as p — > 0, \<j>\ behaves as p k (thus it vanishes on the z— axis, unless 
k = 0). At the same time, a is always nonzero at p = (where \cr\ takes its maximal value) and vanishes as 
p — > oo. (This is why <j> and a are usually refered to as vortex and condensate fields, respectively.) Then the 
phase of cf> changes by 2irk after one revolution around the vortex, while the phase of a increases along the 
vortex. 

Solutions with this behaviour can be constructed by using a standard ordinary differential equations 
solver; typical profiles for \<f>\, \a\ and a discussion of their basic properties can be found in [38]. These 
configurations possess a nonzero momentum along the vortex, with T* = k z w\a\ 2 . Also, their energy density 
is finite everywhere. However, for k ^ 0, the total mass per unit length exhibits a logarithmic divergence 
originating in the 0-sector of the model 14 . 

The heuristic construction of a vorton involves taking a piece of length L of such a vortex and identifying 
its ends to make a loop (see Figure 2). Therefore the central axis of the vortex where the field <f> vanishes 
(for k 0), becomes a circle of radius R — L/2ir. Then the coordinate z along the vortex becomes periodic, 
being replaced by the azimuthal angle ip. This implies that the constant k z in (3.51) also becomes quantized, 
k z = 2nn/L. In such a heuristic construction, the momentum T z along the vortex will then circulate along 
the loop, thus implying a nonzero angular momentum density. Then the corresponding centrifugal force will 
compensate the tension of the loop, thereby producing an equilibrium configuration. 

The vorton would inherit some basic properties of the initial vortex configuration. For example, for 
k 7^ 0, the magnitude of the field (f> would vanish on a circle of radius R in the equatorial plane (the initial 
position of the vortex). At the same time, the maxima of the condensate scalar field a would be located 
on the same circle at 9 = n/2. Also, it is clear that a configuration generated in this way would possess an 
intrinsic toroidal symmetry, the angle ip for the vortex becoming an S 1 generator of the resulting torus, see 
Figure 2. 

14 Notc that this divergence is cured when the <j>— field is gauged. 
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Based on this heuristic construction, a description of large radius vortons has been proposed in [10], 
within an approximation of the vortons' cross section by that of an infinite string carrying a current 15 . 

From the above analysis, it is clear that the resulting toroidal configurations can be described by the 
following axially symmetric ansatz, 

<t> = \4>{r, 9)\e i ^ r ' e) - X(r, 9) + iY(r, 9), a = \a(r, 0)\ e ^+ wt ) = Z(r, 0) e ^+ wt ) . (3.52) 

Then any axially symmetric vorton configuration can be characterized by a pair of integers. Apart from 
the azimuthal winding number n which enters the ansatz (3.52) explicitly, they inherit 16 also the integer k 
of the initial vortex configuration (3.51). This can be seen as follows: the amplitude of the scalar field <f> 
should vanish for some value of r at 6 = n/2, while its phase twists by 2nk after one revolution around the 
boundary of the half-plane consisting of the z— axis (with z = rcos9) and a half-circle at infinity. 

Then, following [14], one can formally define a 'topological charge', which is the product of the two 
integers characterizing the S 1 x S 1 topology of the torus 

N = nk. (3.53) 

However, one should emphasize that, for the model considered in this paper, the solutions do not describe 
topological solitons and thus (3.53) is not a topologically invariant quantity. TV = nk becomes a topological 
charge only in the the sigma-model limit (2.11) (see Ref. [7]). 

Unfortunately, in the absence of exact solutions, it is extremely difficult to follow the details of the 
above heuristic construction 17 . Moreover, while the numerical construction of the initial vortices does not 
pose special problems (since one deals with ordinary differential equations), the systematic study of vortons 
represents a numerical challenge 18 . An approximate form of the solutions can be found on the boundaries 
of the domain of integration only. For example, the leading terms in the asymptotic form of the solutions, 
as r — > oo, read (restoring the dimensionful variables) 

* = ^ + -> H = 1 --^^ 1 + 3cos2e ) + --- + r^ r + ---' H = 7^ r + -' ( 3 - 54 ) 

with X = \<j>\ cos tpi,Y = \<f>\ sin ipi and A, B, C arbitrary constants. Note, that finiteness of the mass imposes 
an upper bound on the frequency, 

w 2 < w 2 max =m%~\ Wa = Ml . (3.55) 
When expressed in terms of the dimensionless parameters, the above condition reads 

ft > w 2 + iftft 2 , (3.56) 

a form which helps to understand the critical behaviour of the numerical solutions. 

Also, as proven in [14], vortons in a flat spacetime background satisfy the virial identity 

T + U=^-Q (3.57) 
T = J d 3 x(\Vcj)\ 2 + |Vcr| 2 ) , U = J d 3 x U(\<p\,\a\), Q = 2w J d 3 x \a\ 2 , (3.58) 



with 
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15 However, as discussed in [18], this approximation fails to provide a good quantitative description, at least for the vorton 
solutions investigated in that work. This is presumably due to the fact that the initial vortices have an infinite energy per unit 
length. However, a similar description of large radius vortons provides excellent agreement with the full numerical simulation 
for d = 2 + 1 dimensions [42] . 

16 Note, that the number k can be defined also for non-axially symmetric configurations. 

17 Remarkably, a similar heuristic construction could be demonstrated in d = 5 Einstein gravity [41]. The large radius black 
rings there are essentially boosted black strings. 

18 The correspondence between vortices and vortons could be studied in detail for planar kinky vortons, i.e. solutions in 2 + 1 
dimensions [42]. The starting configuration represents a kink string, which is known in closed form. This leads to a number of 
exact results and greatly simplifies both the analytical and numerical treatment of the problem. 

19 Here the operator V is expressed with respect to the three dimensional cuclidcan metric. 
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Figure 3. The scalar functions X, Z, the energy density — Tjf and the angular momentum density T,£ are shown for 
a typical k — 'semitopological' vorton. The left plots are for the z — plane, while the right plots are for a plane 
containing the axis of symmetry. The parameters of the solution are w — 0.75, n — 2 and f3\ = 0.71, P2 = 2.22, 03 = 
1.7. 
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Figure 4. Le/t- The mass of k = 'semitopological' vortons in a fiat spacetime background is shown as a function 
of the frequency for several values of the azimuthal winding number n. The vertical line corresponds to the maximal 
value of the frequency as given by (3.55). Right: The M(Q) diagram for the same solutions. Also shown is the mass 
for Q free bosons, M = M a /M<f,Q as given in units set by the scalar field 0. Note, that in all figures in this work 
exhibiting results for families of solutions, the large dots represent the data points. 



which clearly shows that the solutions owe their existence to the harmonic time dependence of the condensate 
field a. Note, that (3.57) holds for generic solutions of the model, without any symmetry assumption. The 
mass of the solutions can be written in terms of the above quantities as 

M = T+-wQ + U. 

For spinning axially symmetric solutions, a more suggestive form of (3.57) is found by replacing Q = J/n. 
Then one can see that the angular momentum of the solutions is bounded from below, i.e. the limit J — > 
is not allowed. The relation (3.57) is also used to verify the accuracy of the numerical results. 

A rather similar virial identity holds for Q— ball solutions [14]. These are solitonic solutions in a model 
with a single complex scalar field which vanishes at infinity (like cr) and possesses a non-renormalizable 
potential. This suggests to view the toroidal solutions as <fi— dressed Q— balls made of the complex er-field, 
where the non-renormalizable |<t| 6 term in the Q-ball potential is replaced by the interaction with the scalar 
field 4>. As we shall see, the numerical results give some support for this analogy. However, the vortons 
possess also some special features which singularize them. 

3.2 k = 'semitopological' vortons 

These are the simplest solutions of the model studied in this work. In a heuristic construction, they are 
found by starting with loops made of vortices possessing k = in the vortex ansatz (3.51). This implies that 
the scalar field is purely real, i.e. Y = 0. 

Moreover, since (f> does not vanish at the core of the vortex, then, for the corresponding axially symmetric 
solution, X should be a nodeless function 20 . 

The azimuthal winding number n is non-trivial in this case, but the 'topological charge' N — kn as 
defined by (3.53) is zero, and so strictly speaking these solutions should not be called vortons. Then in what 
follows we shall call these configurations 'semitopological' vortons. 

Also, in some sense they represent axially symmetric generalizations 21 of the spherically symmetric 

20 Indeed, for all solutions with a real <j> which we could construct, the function X(r, 9) was strictly positive. 
21 The FLS model [22] is recovered for a vanishing cr-potential in (2.2). However, we have found that the basic properties of 
the more general solutions here hold also in that limit. 
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Figure 5. The dependence of the k — solutions on the parameters fa is shown for several values of the 
winding number n. 
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Figure 6. Left: The dependence of the k — solutions on the parameter fa, is shown for several values of the 
winding number n. Here and in Figure 5, the vertical lines correspond to the limiting values of fa as given by (3.56). 
Right: The mass of typical k = vorton solutions is shown as a function of the azimuthal winding number n for fixed 
parameters of the potential and a given frequency. 



solitons studied by Friedberg, Lee and Sirlin in [22]. 

For n 7^ 0, the k = configurations satisfy the same boundary conditions as in (2.34), (2.35), (2.36), 
(2.38) (with Y = this time). The profiles of a typical solution are shown in Figure 3, together with the 
corresponding distributions for the energy density — T t and the angular momentum density T* There we 
plot the quantities for a slice corresponding to the z — plane (left panel) and also for the y — plane (with 
x, y, z the usual cartesian coordinates). 

One can see that the functions X, Z possess a pronounced angular dependence, with a nonzero minimum 
of X located in the ■& = tt/2 plane (z — 0), at some value r = R. There the function Z approaches the 
maximal value, and likewise the energy and angular momentum densities. Based on that, we can define 
formally an effective ring radius R also for the k = solutions 22 . The energy density typically exhibits a 



Alternatively, following [40] one can define a mean radius of the configurations R = , 3 
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torus-like structure. (This becomes apparent by considering surfaces of constant energy density of e.g. half 
the respective maximal value of the energy density.) 

Here we note that the k = configurations are obtained easily, and with a better numerical accuracy 
than any other solutions in this work. Therefore it was possible to investigate some regions of the five- 
dimensional parameter space in a systematic way. Starting with the dependence on the frequency, we note 
that the behaviour found in [16] for flat spacetime spinning Q-ball solutions is recovered, see Figure 4. 
(Similar results have been found also for other sets of input parameters.) As one can see, the solutions 
exist for a limited range of frequency, w m i n < w < w max (with w max = M a jM^, for the scaling used in this 
work, see (3.56)). Both mass and charge (or, equivalently J = nQ) diverge at the limits of this interval. 
They assume minimal values at some intermediate value w cr . Then a diagram M(Q) shows the existence 
of two branches, see Figure 4 (right). For the same Q, the solution with w < w cr has lower mass than the 
corresponding one with w > w cr . Moreover, those configurations which have a smaller mass than a collection 
of non-interacting field quanta with the same charge are expected to be stable with respect to decay to the 
free particles 23 , at least for small values of n. 

The dependence of the solutions on the parameter j3\ is exhibited in Figure 5 (left). One can see 
that the solutions exist only for a limited interval, /Si( m j n ) < /?i < flir m ax)> an d the mass diverges as the 
limits are approached. Interestingly, the value of f3ir ma x) does not depend on the winding number n, with 
Pi( ma x) = a/2(/3 3 - w 2 )//?2, see the relation (3.56). 

When studying instead the dependence on /3 2 , one first notices the existence of axially symmetric gener- 
alizations of the FLS spherically symmetric solitons [22], which are found for /3 2 — > 0. Also, the mass seems 
to diverge when a critical, maximal value of /3 2 is approached, see Figure 5 (right). As implied by (3.56), 
this maximal value is given by 2(/?3 — w 2 )/ '/3 2 (and thus again it does not depend on n). 

We have performed a similar study also for /3 3 , see Figure 6 (left). One can notice the existence of a 
minimal value of this parameter, iiz(min) = w2 + , see (3.56). As ^3 — > f3 3 ^ min ^ the mass diverges. At 

the same time, no upper bound on f3% appears to exist, at least for the set of parameters we have considered. 
Finally, in Figure 6 (right), we plot the dependence of mass on the winding number n (where one can notice 
an almost linear relation). Also, as expected, the angular dependence of the scalar functions increases with 
the winding number n and the location of the biggest angular splitting moves further outward, as well as 
the global maximum of the energy density. 

Note also that for all cases in Figures 5, 6, the charge Q shows a similar pattern. (Here we recall that in 
all plots in this work the global charges are expressed in natural units set by r]^,X^.) 

The k = solutions possess also an interesting spherically symmetric limit, in which case n = and 
X(r, 9) = X(r) and Z(r,&) = Z(r) (with Z(0) 7^ this time). Some properties of these static solutions are 
discussed in a more general context in Appendix C. Here we mention that, as seen in Figures 4-6, they share 
the basic properties of the spinning solutions. 

We close this part by remarking that since the study of the k = solutions does not pose special numerical 
difficulties, and also the corresponding vortices possess a finite mass per unit length, these may open the 
possibility to work out the details of the heuristic construction described in Section 3.1. We hope to report 
elsewhere on that. 

3.3 k — 1 vortons 

The distinguishing feature of these solutions is that the vortex field 4> has a zero on a circle in the 9 = tt/2 
plane (i.e. at z = 0). The profiles of a typical solution are shown in Figure 7. There the plots are again 
given for the (xy) and (xz) planes. (Note, that the X, Y, Z, |</>| profiles are shown in Figure 15 in terms of 
cylindrical coordinates (p, z).) We also recall that x, y, z correspond to dimcnsionless coordinates, the length 
scale being set by l/M^). 

A description of the generic properties of the scalar functions can be made as follows. At r = 0, Z 
vanishes, while X has a local (negative) extremum. Near the origin, the functions X, Y, Z show a large 

23 This conjecture is based also on similarities we have found between the k = solutions and the spinning Q-balls in a model 
with a single complex scalar field. In that case, as discussed in [43] (though for a different formulation of the problem), the 
n = 0, 1 solitons arc stable for some range of frequency, both at the linear and non-linear levels. 



16 




Figure 7. The scalar functions X, Y, Z and the amplitude V X 2 + Y' 2 of the field <j> are shown for a typical k = 1 
vorton solution. 
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Figure 8. The energy density — T* and the angular momentum density are shown for the same solution in Figure 
7. The input parameters of the solution are w = 0.05, n — 2 and /3i = ySa = 1, /fe = 0.51. Here and in Figure 7, the 
left plots are for the z = plane (except for the F-function, where we take z — 80), while the right plots are for a 
plane containing the axis of symmetry. 

variation, with Y approaching its extremum. Then there follows an intermediate region located inside the 
vortex, where X stays very close to zero, while the function Z is almost constant and close to its maximal 
value. (Typically, this region is on the order of 10 2 units of the fundamental length scale 1/M$.) There, the 
function X crosses the value zero on a circle of radius R at z = 0, with a small but non-zero first derivative, 
d r X(R, 7r/2) ^ 0. Finally, there is a third region outside the vortex ring, where the fields approach their 
vacuum values. There, because of (3.54), the function Y has a slower decay than the functions X and Z. 
Also, both X and Z show large variations in a transition domain between the second and the third region, 
approaching afterwards very fast their asymptotic values. 

The corresponding plots for the energy density — T\ and the angular momentum density T* are shown in 
Figure 8. One can see that, as expected, the shape of copies the profile of the scalar function Z. However, 
the energy density distribution has some special features which single out vortons among the field theory 
solitons. For all studied configurations, the energy density has a local maximum at the origin , whose value 
decreases with n. Most of the vorton's energy is concentrated in a shell associated with a circular domain 
wall of cross-section tq. (The value of ro gives the transition between the second and third characteristic 
regions as defined above.) The center of this shell is located at r = R, 6 = ir/2. Then, for the typical 
solutions in this work, the shape of the vorton energy density resembles a hollow tube (although there is 
always an almost constant nonzero value of T\ inside that tube). When the azimuthal winding number n 
increases, the position of this tube moves outwards. More relevant plots for typical solutions are given in 
Ref. [14] (see also the profiles in Figures 15, 16). 

The vorton topology for a solution with k — 1 is illustrated by the diagrams in Figure 9. On the left 
picture we plot the lines of constant amplitude for the scalar <f> (closed curves encircling the center of the 
vortex where <fi vanishes) and phase (radial lines emanating from the center of the vortex) for a typical vorton 



Note that, for spinning Q-balls, the energy density vanishes at r = 0. However, this is not the case for vortons, due to the 
nonzero contribution to Ti of the </<-field potential. 
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Figure 9. The behaviour of the field <f> is displayed for a typical k — 1 vorton solution with /3j = 1, /?2 = 0.97, 
03 — 0.55, to = 0.13 and n — 2. In the left panel we plot levels of constant amplitude \<j)(x, 0, z)\ (closed lines) and 
phase ip(x,0, z) (radial lines). In the right panel, the horizontal and vertical length of each arrow is proportional to 
the real and imaginary part of <j> (X and Y , respectively). The red points there indicate the position of the ring. 

solution. The same picture is shown on the right panel in a different representation. There, the horizontal 
and vertical length of each arrow is proportional to the real and imaginary part of (f> (X and Y, respectively). 
One can see, for instance, that the function Y takes significant values only in a region with a radius of about 
50 /Mtj,, surrounding the vorton's core. 

A systematic analysis of the parameter space of vorton solutions would be a very difficult task that we 
did not aim at here. Also, unfortunately, some regions of the parameter space are not accessible within our 
numerical approach. The solutions investigated so far cover a domain centred around the point 25 {/3\ = 1, 
P2 = 0.95, /?3 = 0.5}, with maximal deviations of /?,; around these values of about 7 percent. Also, all 
configurations reported here (including the gravitating solutions discussed in Section 4.2 below) represent 
thick vortons. Thus the radius of the vortex ring does not differ significantly from the thickness of the vortex. 
The constructions of thin vortons seems to require large values for the winding number n. However, we did 
not yet manage to find within our numerical approach solutions with n > 6 with high numerical accuracy. 

We also mention that these vortons have been constructed by descending from solutions of the sigma 
model limit (2.11), where the constraint X 2 + Z 2 + Y 2 = 1 is satisfied. However, most of the solutions in 
this work exhibit large deviations from this value, with \(j)\ 2 + \a\ 2 ranging typically from about 0.8 to about 
1.1. 

The picture we have found when attempting to vary the input parameters seems to be rather similar to 
the one found in the k — case. (Although, for true vortons, the study of the limiting solutions is much 
more difficult.) In particular, the dependence of the solutions on the frequency looks qualitatively similar 
to the one exhibited in Figure 4. The Q— ball type behaviour is recovered. The solutions exist for a limited 
range of frequencies, and M(Q) exhibits again two branches. (Plots illustrating this behaviour can be found 
inRef. [14].) 

Provided both k = and k = 1 solutions exist for a given set of parameters, then the mass and angular 
momentum of the k = 1 solutions are found to be usually larger than those of the corresponding k = 
configurations. (Note, that we could construct k — 'semitopological' vortons for a much larger set of values 
of the potential parameters Pi than in the k = 1 case.) 

An interesting feature of the k = 1 solutions which, to our knowledge, has not been reported so far in the 
literature is their nonuniqueness. This means that for the same input parameters (f3i,n,w) we have found 
two vorton solutions with different global charges. The situation is illustrated in Figure 10, where we plot a 

25 The thin vorton solutions in [18] have been constructed for a choice of the parameters of the potential fi\ = 1, /?2 = /?3 = 2/3. 
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Figure 10. The profiles of the scalar functions X(r,8), Y(r,0), Z(r,6) and of the energy density — T| are shown 
for two different k = 1 vorton solutions with the same input parameters: /3i = 1, 02 = 0.97, /?3 = 0.54, w = 0.135 
and n = 4. The solution with dotted lines has higher mass and angular momentum and is an excitation of the 
fundamental solution, shown with solid lines. 

fundamental solution, with lower mass (continuous line) together with the corresponding excited one (dotted 
line). The difference between the masses and charges of these two solutions increases with the winding 
number n. For the parameters in Figure 10 one finds e.g. M( fund ^/M^ exclt ^ ~ 0.98 for n = 2, while 
M (fund.)/ M (cxcit.) ^ q 9 for n _ 5 interestingly, the solutions with higher mass have a smaller radius 
that the lower mass (and thus fundamental) configurations. Also, the excited configurations display a more 
complicated behaviour of the function X at 8 = 7r/2, exhibiting always an intermediate local extremum 
there. The profiles of the energy densities for the fundamental and excited solutions look also very different, 
see Figure 10. (Note also, that in Figures 7, 8, 9 we have displayed fundamental solutions; this holds also 
for all configurations reported in [14]). 

These excited states are expected to be unstable. They are much more difficult to construct than the 
fundamental solutions 26 . For the input parameters investigated so far, all fundamental solutions appear to 
possess excited configurations. However, the generic picture is likely to be more complicated, and other ex- 
cited solutions may exist as well. For example, the nonuniqueness of solutions for the same input parameters 
and the same topological sector has been noticed in [44, 45] for the case of Yang-Mills-Higgs theory, with a 
Higgs field in the adjoint representation. As shown there, new branches of solutions appear at critical values 
of the strength of the Higgs self-coupling parameter A. 



26 The excited solutions have been found as limiting configurations of a secondary branch of gravitating solutions, see the 
discussion in Section 4.2. 
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Figure 11. The profiles of the scalar functions X(r, 8), Y(r, 9), Z(r, 6) and of the energy density — T\ are shown for 
a typical di-vorton with j3\ — f3% — 1, f3z = 0.5, w — 0.023, n = 1. 

In this work we have only studied the dependence on the frequency of the excited solutions for fixed /3j,n. 
The results show that the picture found in [14] for the fundamental configurations holds also in this case. 

Also, it is not clear to us how to interpret the nonuniqueness result in terms of the heuristic construction 
discussed in Section 3.1. This would imply the existence of excited vortex states for the same input data. 
(The corresponding ansatz is given by (3.51) with k = 1.) However, so far we could not construct such 
configurations. 

3.4 Composite configurations. The di-vortons 

As noted in the introduction, there is a striking analogy between the heuristic construction of a vorton 
(starting with a vortex) and that of a black ring (starting with a black string). However, in the gravity case, 
a variety of solutions describing composite configurations (e.g. a Saturn [46] and a di-ring [47]) have been 
constructed inspired by rather similar arguments. Thus it is reasonable to expect that the global Witten 
model (and, in general, any field theory model admitting vorton solutions) would also possess composite 
solutions-the analogs of d = 5 gravitating Saturn, di-ring and further multi-component objects. 

In what follows we give some preliminary numerical evidence for the existence of a solution describing 
two concentric vortons, i.e. a di-vorton. In a heuristic construction, this configuration is found by taking 
finite pieces of two parallel vortices and bending them in the same plane, around a common center. Again, 
the resulting configuration would be supported against collapse by the centrifugal force. 

The profile of a typical di-vorton solution is shown in Figure 11. There one can see that \(f>\ vanishes 
in the equatorial plane 8 = 7r/2, for two different values of r. These values correspond to the positions of 
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the two concentric rings. Also, the functions X, Y, Z show a much more complicated dependence on r and 
9 than in the single vorton case. One can see e.g. that, different from the single vorton case, the function 
X starts from a positive value at r — (which is close to one, for all solutions we could construct), takes 
negative values for some range of r and approaches the v.e.v. X = 1 asymptotically. The mass and angular 
momentum of these solutions are much larger than the mass of two (usual) vortons with the same values 
of the input parameters. (Note also, that all di-vorton solutions we could construct so far correspond to 
configurations close to the sigma- model limit.) The angular momentum of the solutions is concentrated in 
two tori centred around the locations of the rings 27 , while the mass density exhibits a local extremum at 
r = 0. The two radii of the di-vorton are also much larger than in the single vorton case. 

The numerical accuracy was much lower in this case, which has prevented us from a systematic study of 
di-vortons. However, we expect them to share the basic properties of the usual vortons. In particular, they 
should exist for a limited range of frequencies. 

4 Gravitating vortons 

When a is increased from zero, while keeping (A;n, w) fixed, a branch of rotating solutions emerges from 
the flat spacetime configurations. Since the dependence on a is different for k = 0, 1 solutions, we shall 
discuss them separately. 

4.1 Gravitating k = 'semitopological' vortons 

Not entirely surprising, some properties of the gravitating k = solutions are similar to those of the even- 
parity rotating boson stars in a model with a single complex scalar field [16]. 

The behaviour of a typical k = solution is exhibited in Figure 12. These two-dimensional plots exhibit 
the small-r dependence for three angles 9 — 0, 7r/4 and ir/2. (We recall that these solutions are symmetric 
w.r.t. an equatorial plane reflection.) One can notice, for example, that the functions /, I and m exhibit a 
strong angular dependence, with I ^ m (except on the z— axis). Also, the metric functions are completely 
regular and the line element shows no sign of a horizon. However, the shape of the scalar functions X, Z is 
similar to that found for solutions in a flat spacetime background. 

Considering the spatial structure of the solutions, we have found that the scalar fields X, Z typically give 
rise to a torus-like energy density, just as for the flat spacetime solutions. (This torus-like shape becomes 
apparent by considering surfaces of constant energy density.) Again, the maximum of the energy density of 
the solutions increases and shifts towards larger values of the coordinate p — rsm9, when n (and thus the 
angular momentum) increases. 

To demonstrate the effects of gravity, we exhibit in Figure 13 the mass and angular momentum of typical 
solutions as a function of frequency. Three different values of a are considered there, for the same parameters 
Pi of the potential and a winding number n = 1. Similar to the case of a flat spacetime background, the 
solutions exist for a limited range of frequencies, w m i n < w < w max , with w max given by (3.55). However, 
as expected, the divergencies of the global charges at the limits of the w— interval which were found for 
a = (see Figure 6) are regularized by the effects of gravity. Moreover, for q/0, the solutions exhibit 
a complicated pattern after approaching w m i n . There one observes an inspiralling behaviour of the k = 
solutions, towards a limiting solution at the center of the spiral, for a frequency wu m > w m i n . (This 
behaviour could be studied more precisely for spherically symmetric solutions (n = 0), which we do not 
display here.) As a result, when the mass M is considered as a function of the scalar charge Q, we observe 
a cusp structure, as illustrated in Figure 13 (right). This is essentially the behaviour noticed in [16] for 
rotating boson stars. The branches starting from the maximal value of the frequency and ending at the 
maximal values of the mass are expected to be stable, in analogy to boson stars. 

The dependence on the parameter a is shown in Figure 14. (A similar behaviour has been found for 
spherically symmetric solutions, n = 0.) One can see that the solutions exist for arbitrarily large values of 



27 Different from the case of a Saturn [46] or of a di-ring [47] in d = 5 Einstein gravity, both components of a di-vorton carry 
an angular momentum density with the same sign, see the fiat spacetime limit of the relation (2.31). 
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Figure 13. Left: The mass of gravitating k = 'semitopological' vortons is shown as a function of the frequency 
for several values of the coupling constant a. The vertical line corresponds to the maximal value of the frequency as 
given by (3.55). Other input parameters of these configurations are /3i = 0.03, /?2 = 1, A = 2.5, n = 1. Right: The 
mass-charge diagram for the same solutions. Also we show the mass for Q free bosons, M = M a jM^Q in the units 
employed here. 




Figure 14. The mass and angular momentum are shown as functions of a for gravitating k = 'semitopologi- 
cal' vortons. 

a (i.e. an arbitrarily large v.e.v. of the scalar field <fi). 

This strongly contrasts with the picture found for other models featuring scalar fields with a symmetry 
breaking potential, e.g. monopoles or sphalerons, in which case a maximal value of a is found (see the 
review [21]). In the large a limit, the function Z vanishes, while the function X approaches a constant value 
everywhere, X — > 1, corresponding to the v.e.v. of the scalar field. 

To understand this behaviour, we consider the scaled functions 

X= — — , and Z=-, (4.59) 
a a 

Then, after making this replacement in the field equations and taking the limit a — ¥ oo, the functions X, Z 
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Figure 15. The scalar functions X, Y, Z, the amplitude V X 2 + Y 2 of the scalar field <j> and the metric functions 
/, I, m, W axe shown for a k = f gravitating vorton solution with /3i = f, /3a = 0.97, /?3 = 0.54, w — 0.11, n = 3 and 
a = 0.2. The coordinates here are p = rsin# and z = rcos9. 
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Figure 16. The energy density — T* and the angular momentum density T£ are shown for the same solution 
as in Figure 15. 



3000 



4000 
















\ 3000 










n=3 \ 




















% \ 2000 












0.0675 


0.075 


























n=2 


**»^ 




n=4 











4000 


■ 3600 








3500 




n= t» 


3000 


■ 3400 






S 


28000 30000 32000 




2000 




St'*"*!***" 


n=2 


1000 




A' 


P 1= l p 2 =0.97 p 3 =0.55 a=0.2 



0.12 0.14 
W 



0.18 



5000 10000 15000 20000 25000 30000 35000 



Figure 17. Left: The mass of gravitating k — 1 vortons is shown as a function of the frequency for several 
values of the winding number n. Other input parameters of these configurations are /3i = 1, (3% = 1, /J3 = 0.54 and 
a = 0.2. Right: The mass-charge diagram for the same solutions. 



decouple, the scaled matter fields being described by the Lagrangian 



L s = d^Xd^X + d^*d»a + M 2 a \a\\ 



(4.60) 



(with a = Ze l ( nip+wt }). It is clear that X — is the only physical solution, and we recover a boson star 
model for the complex scalar field er possessing a quadratic potential. A similar behaviour has been found in 
[16] for gravitating Q-ball solutions with higher order terms in the potential: for large values of the coupling 
constant a, all higher order terms in the scalar field potential become irrelevant. 



4.2 Gravitating k = 1 vortons 

Again, these solutions smoothly emerge from configurations in a fixed Minkowski background. The backre- 
action is included by slowly increasing the value of a. 

The 3D profiles of a typical solution are shown in Figure 15 as functions of the (scaled-)coordinates 
p = r sin 9 and z = r cos 9. One can notice the nontrivial behaviour of the metric functions, which strongly 
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Figure 18. The mass M, radius of the ring R and the values at the origin of the metric function m and of 
the scalar function X are shown as functions of the parameter a for gravitating vorton solutions with /3i = 1, /?2 = 
0.97, fis — 0.54, w — 0.13 and several values of the winding number n. One can notice the existence of two branches 
of solutions connecting different flat spacetime (a = 0) limiting configurations. 



deviate from the flat spacetime values. The functions /, I and m approach their minimal value in the equa- 
torial plane at the location of the ring where X = 0. For all solutions studied so far, we did not notice a 
special behaviour for the functions X, Y, Z in the gravitating case as compared to the flat spacetime limit. 

The same holds for the distribution of the energy and angular momentum densities. In particular, as 
seen in Figure 16, the vorton energy density still resembles a hollow tube. (Note, that the local peak at 
r = of the energy density is much smaller for this n = 3 solution than for the configuration with n = 2 
exhibited in Figure 8.) 

As expected, new effects occur when considering the vortons' dependence on the frequency. This is 
shown in Figure 17 for a typical family of gravitating solutions emerging from a flat spacetime fundamental 
configuration. One can see that some features found for k = 'semitopological' solutions are likely to hold 
also in this case. For example, the range of w is bounded, both M and J approaching a maximum close to 
a minimal value of the frequency w m i n . Furthermore, as for the k = solutions, we notice the existence at 
Wmin of a backbending towards larger values of w, leading apparently also to an inspiralling of the solutions 
towards a limiting configuration. Unfortunately, despite our efforts, numerical accuracy currently does not 
allow us to continue the construction of the (expected) spiral in a reliable way. 

This holds also for the issue of a maximal frequency, in which case we could not construct the counterparts 
of the curves in Figure 13, with a fast decreasing mass and charge. The situation here is puzzling, since, 
for all considered sets of input parameters (ft; n; a), we could not construct k — 1 gravitating vortons with 
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frequencies below a critical value w c . As w — > w c , the solver diverges, or, more often, provides a solution of 
the problem with Y — 0. At the same time, we do not notice any special features of the k = 1 configurations 
close to w c ; in particular the metric functions do not exhibit any pathological behaviour. Interestingly, the 
k = configuration which is found in this way possesses almost the same mass and angular momentum as 
the k = 1 vorton solution with w — > w c (although other relevant parameters (e.g. the value of the metric 
function / at the origin) are different). Then one may speculate that the k = 1 branch would continue for 
w > w c with a branch of k = configurations towards w max given by (3.55), where the vacuum state is 
approached 28 . However, due to limitations, resulting from our numerical treatment of the problem, we could 
not clarify this point. 

A distinct feature of the k = 1 gravitating vortons is their dependence on a. The emerging picture 
which results from the study of several different sets of solutions can be summarized as follows 29 : given 
a fundamental flat spacetime configuration with a set of input parameters (/3,;n, «;), the curved spacetime 
configurations are found by increasing the parameter a. The emerging branch of solutions exists up to a 
maximal value of a (which depends on the input parameters). 

As a — > a max all relevant quantities stay finite. (In particular, we do not see any sign for the emergence of 
an event horizon.) At a max a second branch of solutions is found, extending backwards in a. Asq-> 0, this 
branch approaches the excited flat spacetime vorton discussed in Section 3.3. These features are illustrated 
in Figure 18. (We mention that the ring radius R, the absolute value of the scalar function X and the value 
of the metric function m at the origin are smaller for solutions on the excited branch.) 

The fact that there is an upper bound on a for k — 1 solutions can be understood as follows. A similar 
scaling to (4.59) (this time also with Y = Y/a) shows that the limit a — > oo would imply necessarily X = 
(thus X = 1), Y = 0. However, these values cannot be continuously approached, since \4>\ should always 
vanish on a circle in the equatorial plane. In other words, since k = and k = 1 correspond to different 
'topological' sectors of the theory, no continuous transition between them is possible 30 . 

As a result, the solutions exhibit instead a maximal value of the parameter a, with a backbending towards 
a second a = solution in a flat spacetime background 31 . 

Thus we have reached the important conclusion 32 that the effects of gravity impose an upper bound on 
the v.e.v. r]^ of the scalar <j> for solutions to exist. This upper bound depends on the parameters of the 
model, with a general expression 

with c a parameter of order one, whose precise value depends on the other input data. 

Rotating objects in general relativity may possess ergoregions. However, unlike the case of black holes, 
the ergoregions of regular objects like those in this work would signal the presence of an instability [49]. 

Within the metric ansatz (2.21), the ergosurfacc is defined by the condition 

9tt = -f + jW 2 sin 2 6 = 0. (4.62) 



28 Note also that, since \<f>\ vanishes on a circle in the equatorial plane, the vortons cannot approach in a continuous way the 
vacuum state with X = 1, Y = 0. 

29 However, we do not exclude that a different picture may exist for other sets of parameters, not considered in this work. 

30 It may be interesting to note, that wc could not construct solutions in which the function Y takes very small values and 
can be considered as a perturbation around a configuration featuring only the X, Z functions. Related to that, for all solutions 
with 7^0, the function X possesses a node in the equatorial plane. 

31 We would like to mention that the physical meaning of the second limiting solution is not obvious apriori. Since a — > can 
be reached also when taking rj^ — ¥ (see (2.47)), the limiting solution might represent a scaled boson star in a model featuring 
only the scalar field a coupled to gravity. However, as seen in Figure 18, this is not the case, at least for the solutions in this 
work. 

32 One should also mention that, different from most of the other gravitating solutions with a scalar field featuring a v.e.v., 
we did not find any indication for the existence of black hole solutions. (Typically, they emerge from the globally regular 
configurations for a critical value of a [21].) This property should be attributed to the Q-ball features of the vortons induced 
by the field a. In fact, for the simplest case k = 0, n = 0, one can prove the absence of black hole solutions following the 
arguments in [48] . A similar result may hold also in the spinning case. Note, that this is not obvious, since the approach used 
in [48] cannot be easily generalized to the rotating case. In fact, spinning black holes with complex scalar fields posscsing a 
harmonic time dependence have been constructed recently in [32] for a five dimensional anti-de Sitter spacetime background. 
(They require a fine tuning between the frequency of the field and the event horizon velocity.) 
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Figure 19. The location and shape of the vortons' ergoregion is exhibited for solutions with three different frequencies 
(left) and three different values of a (right). The parameters of the potential are j3\ = 0.14, /?2 = 0.974, /?3 = 0.54 
while a = 0.2, n = 3 (left) and w = 0.07, n = 2 (right). 

Then the ergoregion resides inside this ergosurface. Examining the condition (4.62), we find that this 
is not satisfied for most of the gravitating vorton solutions. For example, when varying the frequency, the 
ergoregion typically occurs for a set of solutions close to the backbending point at w m i n , in particular, for 
secondary branch configurations. The solutions with large enough frequencies (in particular those close 
enough to w max ) have gu < 0. (This is again similar to the situation encountered for boson stars [16].) 

When studying instead the influence of the parameter a, we find that an ergoregion may occur for a set 
of solutions close to a max (and with values of w not too far from w m i n ). That is, for the parameter choice in 
this work, the solutions with small values of a (which are expected to be relevant in a realistic context) do not 
exhibit an ergoregion. Some results displaying these features are shown in Figure 19. One can see that the 
vortons' ergoregions also possess a toroidal shape, whose size and location depend on the input parameters. 
Without entering into details, we mention the existence of similar results for k = 'semitopological' vortons, 
where some of these solutions also possess an ergoregion. 

We close this part by remarking that the di-vortons should also possess gravitating generalizations. 
However, due to the severe numerical difficulties, we could construct such solutions for very small values of 
a only. Their systematic study might possibly require a different numerical approach. 

5 Discussion and conclusions 

More than 20 years after vortons have been first proposed, they continue to provide an interesting field of 
research. Posing a severe numerical challenge, the systematic quantitative study of such configurations has 
started only recently [14, 18], within the global version of Witten's U{1) x U(l) theory. 

The first purpose of this work was to give an overview of the properties of vortons in a fiat spacetime 
background. The heuristic construction of these objects as reviewed in Section 3.1 suggests to view the usual 
vortons as the second member (k = 1) of a family of solutions indexed by the number of nodes k of the 
vortex field </>. Then the case k = would correspond to 'semitopological' vortons. These solutions could be 
studied to a greater extent and seem to exhibit most of the basic properties of the usual vortons (k = 1). 
A new result reported in Section 3.3 is the nonuniqueness of the vortons for the same input parameters. 
There we have given numerical evidence for the existence of excited solutions, which play in important role 
when including the gravity effects. Another new type of configuration describing two concentric rings, the 
di-vorton, has been reported in Section 3.4. Here we remark that one expects that Witten's model possesses 
a variety of other composite configurations. For example, it would be interesting to find a solution describing 
a 'Saturn' i.e. a configuration with a central concentration of the energy density surrounded by a ring (this 
would correspond to different zeros of the scalar field <j>). 
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Another aim of this work was to present a first discussion of how gravity affects the properties of k = 
'semitopological' vortons and k = 1 'true' vortons. Our results can be summarized as follows. First, we 
have shown that, when varying the frequency, gravity imposes an upper bound on the mass and angular 
momentum of these objects. The coupling to gravity gives rise to a spiral-like frequency dependence of the 
mass and charge. Second, we have noticed a different behaviour of the k — and k = 1 solutions as functions 
of the parameter a, describing the coupling to gravity. The k — solutions are essentially 'dressed' Q-balls 
in a more complicated version of the FLS model [22]. (This holds also in the flat spacetime limit.) Different 
from true vortons, the limit a — > oo is allowed in this case. For k = 1, one finds instead two branches of 
solutions joining at a maximal value of a. Finally, we have noticed that, similar to the boson star case, 
vortons may possess an ergoregion, with the associated instability. Here one should say that a systematic 
study of the dependence of the properties of the vortons on the various input parameters is still missing and 
continues to present a numerical challenge. 

As avenues for future research, perhaps the most compelling task would be the study of vortons in 
models including gauge fields as well. The simplest case here would be again the local version of Witten's 
U(l) x U(l ) theory The corresponding gauged vortons could be constructed by using a similar approach to 
that described in this work and will be reported elsewhere. 

However, vortons should also exist in other more complex models, e.g. possessing non-Abelian fields apart 
from scalars. An interesting question to ask in this context is the possible existence of 'springs' [51], i.e. 
vortex loops stabilized by magnetic fields instead of rotation 33 . While 'spring' solutions are unlikely to exist 
in Witten's U(l) x U(l) theory, such configurations have been found in other models featuring non-Abelian 
fields, see e.g. [52, 53]. 

Then it would be interesting to reconsider some of the known toroidal solitons based on the heuristic 
construction of a vorton/ 'spring' starting with pieces of strings. For example, vortex ring solutions are known 
to exist in the Yang-Mills- Higgs model with the Higgs field in the adjoint representation [52]. However, these 
configurations seem to correspond 34 to loops made of Yang-Mills- Higgs vortices discussed in [54] . Finally, 
the possible existence of vortons and 'springs' in the electroweak sector of the Standard Model is perhaps 
the most exciting open problem (see [53, 55, 56, 57] for some work in this direction). 
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A System of differential equations 

The metric functions f,l,m and W which enter the Ansatz (2.21) satisfy the following set of coupled non- 
linear partial differential equations: 



33 Note, that the static dipole black rings can be balanced by immersing them in a background magnetic field [50]. 

34 We thank F. Navarro-Lerida and T. Tchrakian for discussions on that and for sharing with us their unpublished results. 
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(A.2) 
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They are equivalent to the following combinations of the Einstein equations: E\ — E r r —E B e — E£ + ^r-E^ = 
(for (A.l)), E r r + E$ = (for (A.3)), E$ - = (for (A.4)), and £* = (for (A.5)), multiplied with 

suitable factors (e.g. —m for (A.l)). 

The equations for the scalar functions result directly from (2.9) and read: 

V 2 A + 1(V0 • (VI) + yA Qa (, 2 - A 2 - r 2 ) - 7 Z 2 ) = 0, 

v 2 y + I ( vo • (vr) + Qa^ 2 - a 2 - r 2 ) - 7 z 2 ) = o, (A.6) 
v 2 z + i ( vo • ( vz) + ^z (i M „ 2 - z 2 ) - 7( a 2 + y 2 ) - + <?L±^ = . 

In the above relations we have defined 

w2 A 1 d < 2 dA , 1 9 . . flA, 

r 2 Or Or r z sm 6* oO oO 

Also, C/(A, y, Z) is the potential of the scalar fields, which results directly from (2.2) after substituting the 
scalar field ansatz (2.22) 

U(X, Y, Z) = Ja (A 2 + Y 2 - r? 2 ) 2 + ^A CT Z 2 (Z 2 - 2r? 2 ) + 7 (A 2 + Y 2 )Z 2 . (A.S) 
These equations are supplemented with the following constraints 

7 (# ^) ^ - - ^ ((W) a - ^(*0 a ) + js (W) 2 ^(*/) a ) (A.9) 

lsin 2 9 /,„ Tir , 2 1 TJ ,, 2 \ 1 /. , 2cot6> \ 2 /„ cot6> \ TYZsin 2 <9. 
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. d r fdef d r ld e l Ism 1 .Am %m , .... 

, d r gl H j5 —5 -f—drWdeW - — — (d r md g l + d r ldgm) - cot 6 (A.10) 

I ] l 2V f z 2lm m rm 

_ e + 9 df)W + 32jtG ( dr xd e X + d r Yd e Y + d r Zd e Z) = 



which result from the Einstein equations E T r — Eg = for (A.9), and E e r = for (A. 10)). 



31 



B The nonvanishing components of the energy momentum tensor 

The nonvanishing components of the energy-momentum tensor of gravitating vortons are found by substi- 
tuting the ansatz (2.22) in the general expression (2.7) and read 

T r r = ( W + ^f^-U { X,Y,Z)-'^Zi 
r J r z sin I 



m 

2/ 



1 



{d r Xf + (d r Y) 2 + (8 r Z) 2 - - ((OgX) 2 + (dgY) 2 + (dgZ) 2 ) 



{d r XdeX + d r Yd Y + d r ZdgZ) , 



T$ = ( W + ^f^-U(X,Y,Z)-'^ 
r J r 2 sin I 
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and finally, the most important components are 



~T\ = I ((VX) 2 + (VY) 2 + (VZf) + U(X, Y, Z) + —1^-^J-Z 2 + (w 2 - ^W 2 )^r, 
m r 2 sin 9 1 r z } 

_ „ , nW . Z 2 

r J 

corresponding to the energy density and the angular momentum density, respectively. 



(B.2) 



C Spherically symmetric solutions 

These solutions occur as a limiting case of the general configurations. Their line element is found from (2.21) 
for / = 772 and W — 0, and reads 



ds 2 = j (dr 2 + r 2 (d6 2 + sin 2 9d^ 2 )) - fdt 2 , 



(C.l) 



the metric functions / and m depending only on the radial coordinate r. The scalar field ansatz results from 
(2.22) for 77 = and reads 



<f> = X(r), a = Z(r)e iwt . 
The equations greatly simplify in this case. From (A. 6), (A.l), (A. 4) one finds 
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w 2 m 



m 



A CT 777 



Z" + (tt- + -)Z' + -n-Z - j-X 2 Z + ^--Z( v 2 a - Z 2 ) = 0, 



"2t77 r 



/ 



/ 



2 / 



for the scalar fields, and 



/" + + \ ~ j)f + l^Gm [- 2 ^f- + U(X, Z) ] = 0, 

777" + - - — + + lG*Gm ( X' 2 + Z> 2 - ^ + Z) ) = 0. 

T 777 2f 2 \ f 2 f 



(C.4) 
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for the metric functions (with U(X,Z) given by (A. 8) with Y = 0), together with the constraint (which 
results from (A. 9)) 

— + ^T- — -Pl+ 32^G(A' 2 + Z' 2 ) = 0. (C.5) 
m j z rm 2m z 

The approximate form of the solutions as r — > reads (note that Z does not vanish there) 

fir) = fo + / 2 r 2 + 0{r A ), m(r) - m + m 2 r 2 + 0(r 4 ) , (C.6) 
X(r) = x +x 2 r 2 + 0(r 4 ), Z(r) = z + z 2 r 2 + 0(r 4 ) . (C.7) 



All coefficients in this expansion arc fixed by the parameters fo, m , x and z (one finds e.g. x 2 = 
W - 4) + l4) and * 2 = ^ (-w 2 + f (jx 2 - (4f (z 2 - ry 2 ))). The large-r expression of the 

scalar functions X, Z results from (3.54) by taking A = 0. The leading order terms in the corresponding 
expansion of the metric function are 

/(r) = l + A + ^ + ... | m(r) = l-^ + ..., (C.8) 

with /i a parameter fixing the ADM mass of the solutions. 

Some basic properties of these configurations are rather similar to those found for spinning k — solutions. 
In particular, we have found the same qualitative dependence on the parameters w, a as in that case. 
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